High-Fidelity Control in High-Q Resonators 
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We describe how the transient behavior of a tuned and matched resonator circuit may be in- 
tegrated into an optimal control theory (OCT) pulse-design algorithm to derive control sequences 
with limited ringdown that perform a target quantum operation with high fidelity in the presence 
of resonator distortions of the ideal waveform. Extending the control system model to include a 
full description of the resonator transfer function allows the use of high quality factor (high-Q) res- 
onators to increase signal-to-noise ratio (SNR) of detection and overall sensitivity without sacrificing 
the ability to perform a universal set of quantum operations. We provide experimental verification 
by implementing a robust bandwidth-limited OCT pulse on an inhomogeneously broadened solid- 
state organic radical spin system. We also generate example pulses that allow universal control of 
anisotropic-hyperfine coupled electron-nuclear spin systems via electron-only modulation even when 
the bandwidth of the resonator is significantly smaller than the hyperfine coupling strength. These 
results demonstrate how limitations imposed by linear response theory may be vastly exceeded when 
using a sufficiently accurate system model to optimize pulses of high complexity. 
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I. INTRODUCTION 

Inductive coupling to a tuned and matched resonator 
circuit is a standard technique for the manipulation and 
measurement of a quantum system. In a common appli- 
cation, a resonator converts a pulsed electrical signal into 
a pulsed magnetic field and the resulting magnetic field 
produced by the quantum system back into an electrical 
signal. The same principles apply as well to capacitive 
coupling to electric fields. The conversion efficiency is 
determined by a number of factors, including the quality 
factor (Q) of the resonator, defined as the ratio of stored 
to dissipated energy in the resonant circuit. All other 
factors being equal, both the magnetic field strength re- 
sulting from the application of an electrical signal of a 
certain power and the signal-to-noise ratio (SNR) of an 
electrical signal induced by the magnetic response of the 
sample scale as y[Q 

An important consequence of using tuned circuits is 
that any resonator has an associated ringdown time dur- 
ing which energy stored in the circuit dissipates. The 
signal to be observed from the quantum system is typi- 
cally many orders of magnitude smaller than the control 
amplitudes, normally requiring a spectrometer deadtime 
of at least five times the ringdown time before receiver 
circuitry, such as low-power high-gain preamplifiers, can 
be switched on. When moving to a high-Q resonator 
to enable faster control and greater SNR, the relative 
amount of stored energy in the circuit increases, lead- 
ing to longer ringdown and a deadtime that can exceed 
the characteristic phase coherence time (Tg) of the quan- 
tum signal, significantly lowering measurement SNR and 
possibly preventing observation of the quantum signal at 
all. A dispersion of interaction strengths between quan- 
tum degrees of freedom in the sample may also lead to a 
spectral breadth which exceeds the resonator bandwidth 



of frequencies which are efficiently converted, prevent- 
ing the rotation or detection of the entire spectrum with 
a single unoptimized pulse For these reasons it is 
common practice, especially in the field of pulsed elec- 
tron spin resonance (ESR), to intentionally spoil the Q 
to reduce it to a value where these effects are no longer 
significant 

Another challenge inherent to the use of resonant cir- 
cuits for pulse transmission is the distortions they pro- 
duce of an ideal pulse waveform, once again due to the 
finite response time of the reactive elements used to con- 
struct a tuned circuit. For low values of Q (roughly 10 
- 200), the effects of waveform distortions are relatively 
minor, resulting in a typical drop in pulse fidelity of a 
few percent [7(. However, as the Q is increased for ap- 
plications that require high sensitivity the resonator dis- 
tortions of the ideal waveform lead to a pulse action on 
the quantum system that differs drastically from a square 
pulse approximation. 

There are many applications of pulsed induction tech- 
niques that require both high sensitivity and the ability 
to accurately control the state of a quantum system. For 
example, thin film samples have a number of well-known 
applications in chemistry and materials science and have 
drawn recent interest in the field of quantum information 
processing (QIP), where a highly structured 2D geome- 
try maps naturally to node-based proposals for quantum 
computing Studies of these samples would ben- 

efit from high-Q, large filling factor resonators [III, flij . 
Additionally, advances in inductive imaging techniques 
have allowed for sub-micron resolution, pushing the lim- 
its of sensitivity using conventional resonators 
As more progress is made on the development of novel 
quantum devices the need for high sensitivity measure- 
ment will increase, requiring precise control with high-Q 
resonators. 
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Recent work proved that integration of a resonator im- 
pulse response function into an optimal control theory 
(OCT) pulse-design algorithm allows broadband excita- 
tion in pulsed ESR with limited power [l8j . Here, we fo- 
cus on applications that require a universal rotation that 
operates identically on all input states, particularly those 
that require a high Q resonator. Specifically, we describe 
how resonator distortions of an ideal waveform may be 
integrated into an optimal control theory (OCT) pulse- 
design algorithm to derive control sequences with limited 
ringdown that perform a target quantum operation with 
high fidelity. We experimentally verified the method by 
optimizing and implementing a bandwidth-limited OCT 
pulse, robust to static and microwave field inhomogene- 
ity, which allows the observation of the free-induction de- 
cay (FID) of irradiated fused-quartz [2(| in a high-Q 
(w 10,000) X-band (w 10 GHz) rectangular cavity. We 
also consider a specific application of achieving univer- 
sal control of a anisotropic-hyperfine coupled solid-state 
electron-nuclear spin system via electron-only modula- 
tion even when the resonator bandwidth is significantly 
less than the hyperfine coupling strength. Through the 
derivation of example pulses, we demonstrate how the 
limitations imposed by linear response theory may be 
vastly exceeded when using a sufficiently accurate sys- 
tem model to optimize pulses of high complexity. 



II. RESONATOR MODEL AND RINGDOWN 
SUPPRESSION 

Several models of the transient behavior of a resonator 
have been suggested and analyzed in detail [2lTl23j , with 
various methods suggested for the suppression of pulse 
transients [24Tj27| . We use the model of Barbara, et al. 
[22} consisting of a series RLC resonant circuit capaci- 
tively coupled to a voltage source through a real 50 Ohm 
impedance (Figure [I}. We restate here the results of 
their analysis relevant to our integration of the model 
into an OCT algorithm, with further details available in 
the references. 

The reactive elements necessary for efficient signal con- 
version have a finite response time to variations of the 
control signal, causing distortions of the ideal waveform. 
For a linear, time-invariant system, the waveform inci- 
dent upon the quantum system, y(t), is given by the con- 
volution of the ideal control pulse, x(t), with the impulse 
response, h(t), of the transmitting circuit, 

/oo 
h(t - t)x(t) dr. (1) 
-00 

The impulse response gives a complete description of the 
transient behavior of the transmitting circuit and may 
be either inferred from measurement, or derived from an 
appropriate model. The ideal pulse is described as a time 
varying voltage of amplitude v(t) and phase 4>{t) applied 
at a driving frequency, uj t : 



When the driving frequency is set near resonance with 
energy level splittings of the quantum system, the mag- 
netic field induces transitions between the eigenstates of 
the system Hamiltonian. This model is representative of 
a general quantum control scenario and allows us to ana- 
lytically understand the dynamics of high-Q resonators. 

The resonator circuit response to an impulse of con- 
stant amplitude and phase may be determined by the ap- 



plication of Laplace transformation techniques [28j . The 
resulting filtered output waveform in the time domain is 
determined by the inverse Laplace transform, calculated 
via a partial fraction expansion, of the filtered waveform 
in the s-domain. The transient coil current, iz,(f), which 
is proportional to the magnetic field applied to the sam- 
ple, may be represented as a sum of poles, 5k, and corre- 
sponding residues, dk, after inverse Laplace transforma- 
tion: 



i L (t)ocY, d ke~ 5kt - 

k 



(3) 



The four reactive elements in the circuit lead to four 
poles of the transfer function: 

1. A steady-state oscillation at the driving frequency, 

2. An exponentially decaying DC transient. 

3. An exponentially decaying oscillation at the probe 
free-ringing frequency, fi. 

4. An exponentially decaying oscillation at minus the 
probe free-ringing frequency, -Q. 

It is the third and fourth poles that are of most interest 
for determining the transient response of the circuit, as 
the DC transient is normally very small and the steady- 
state oscillation may be eliminated by moving into a ref- 
erence frame rotating at the driving frequency. In a rotat- 
ing wave approximation, this frame rotation also allows 
us to neglect the effect of the fourth pole. 

The relevant third pole, S3, may be written as 



SI = 7 - iuof 1 



4i?o/ 2 



(4) 



where 7 = ujq/Q is the rate of transient decay in terms 
of the inductor quality factor, Q = uoL/r, and 



(5) 



is a scaling factor determining how close the resonator 
fr ee-ringin g frequency, ft, is to the tuned frequency, ujq = 
y/l/LCx- For high-Q resonators, r << Rq, Q ~ u)q, 
and the transient oscillations in the rotating frame are 
critically damped. The resulting dynamics may then be 
approximated as a simple exponential ringup and ring- 
down of the pulse amplitude with time-constant 



V s (t) = v(t) cos (oj t t + 0(f)). 



(2) 



T r = Q/UJQ. 



(6) 
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This time-constant differs from the commonly used 
2Q/ujo, which is valid only for an isolated series RLC cir- 
cuit without the inclusion of a matching capacitor (2lj . 
Such approximations are used here only for the sake of 
argument and demonstration; when high-fidelity control 
of a particular system is desired the full form of the res- 
onator transfer function must be experimentally deter- 
mined. 

Resonator ringdown may be suppressed by the appli- 
cation of a trailing compensation pulse of appropriate 
length and amplitude to drive the energy stored in the 
resonator to zero at the end of the pulse [29|, |30| . As 
shown in Figure [1] this compensation pulse can signifi- 
cantly shorten the spectrometer deadtime, but at the ex- 
pense of introducing an additional rotation to the quan- 
tum system. In the next section, we describe how the 
resonator transfer function and a ringdown compensation 
pulse may be integrated into an optimal control theory 
(OCT) algorithm to enable high-fidelity control even in 
the presence of transient effects and the additional rota- 
tion introduced by ringdown suppression. 



III. OPTIMAL CONTROL THEORY FOR 
HIGH-Q RESONATORS 

In the absence of dissipative processes, the problem of 
control sequence design may be stated as optimizing a 
time-dependent Hamiltonian, H(t), such that the action 
of the Hamiltonian on a general quantum state, p, yields 
the desired dynamics, as governed by the Liouville-von 
Neumann equation, 



dp 
~dt 



lP,H(t)} 



(7) 



As is common practice for dynamical calculations, we will 
use the convention h = 1 for the remainder of this work. 
The formal solution of this equation may be written as 



p(t) = U(t)p{Q)U\t), 



(8) 



where U(t) is a unitary propagator representing time evo- 
lution under the applied Hamiltonian: 



U{t)=Te- l i» H ^ ds 1 



(9) 



where T is the Dyson time-ordering operator accounting 
for non-commumitivity of H{t) with time. For this work, 
we are interested in finding a time-dependent Hamilto- 
nian that accurately implements a desired unitary opera- 
tion, Ud . The notion of accuracy for unitary operations is 
quantified by defining an appropriate performance func- 
tional, such as the average gate fidelity for Hilbert space 
dimension D: 



$(U d ,U pnlse ) = |(C/ d |£/ pulso )| 2 = 



Tr 



t/pulse 



■ 2 D 



(10) 



This performance functional has been shown to be equiv- 
alent to the average state fidelity over a complete set of 



input states 3l| and guarantees the desired evolution 
over all possible input states when $ is close to unity. 

The set of unitary operations that may be generated 
depends on the form of the accessible physical interac- 
tions present in a given system. These interactions are 
commonly separated into drift Hamiltonian terms, Hd, 
which are time-independent, and control Hamiltonian 
terms, H c (t), which may be varied during the course of 
an experiment. We parameterize the control Hamiltoni- 
ans into a set of constant Hamiltonians, {Hk}, for which 
we may independently vary amplitudes, Ufc(i), to yield a 
total system Hamiltonian of 



H{t) = H d + J2Mt)H k . 



(11) 



To avoid having to evaluate a set of time-ordered expo- 
nentials to solve , it is common practice to discretize 
the time-dependence of the control amplitudes into N 
intervals of length tj, such that the unitary propagator 

for a time T = Ylj=i tj ma y ^ e formally expressed as a 
product of time- independent exponentials: 



N 



U(T)=l[U J = l[e 



-i(H d +T, k ulH k ) tj 



(12) 



3=1 



For the sake of simplicity, we will assume tj — At for 
all j. The validity of this time-slicing method of prop- 
agator evaluation depends on the chosen value of At, 
with accuracy increasing for smaller time discretization 
intervals. With all time-dependence effectively removed 
from the problem, the control objective now simplifies 
to determining the set of constant amplitudes, {wj^jL-^, 
which maximize the value of the average gate fidelity 
(fTU|) . This objective may be achieved through applica- 
tion of optimal control theory (OCT) [32|]. In the fol- 
lowing, we utilize the efficient GRadient Ascent Pulse 
Engineering (GRAPE) OCT algorithm [H, which has 
been used extensively in spin-based implementations of 
quantum information processing (see [34 - l36l | for exam- 
ple applications of the GRAPE algorithm) . This method 
is limited to finding solutions that correspond to local 
optima, but high-fidelity pulses may still be consistently 
found, insensitive to the initial guess, that are useful in 
practice [34| . 



A. Optimizing Infinite-Bandwidth Controls 

The relative efficiency of the GRAPE algorithm is de- 
rived from the need to only compute the propagator cor- 
responding to each time step, Uj, a single time at each 
iteration of the optimization. The resulting single itera- 
tion propagators are then used to calculate both the value 
of the performance functional and an approximation of 
gradients used for pulse updating. The truncation of the 
gradients is valid when At is chosen to be significantly 
smaller than the inverse magnitude of the system Hamil- 
tonian. Approximating the gradients in this way avoids 
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the use of finite-differencing and other inefficient meth- 
ods of computing the gradients. In the notation of [33| . 
the average gate fidelity and corresponding gradients to 
first order in At are 



$ = |(E/d|[/ P ui S c(T))| 2 = (PjIXj) (XjlPj) (13) 



and 



<5$ 

— = -2Re [{PjliAtHkXj) (X^)} , (14) 
«K 

where Pj = U] +1 ...U\,Ud and Xj = U r ..U 1 are the back- 
ward and forward representations, respectively, of the to- 
tal pulse propagator at the j th time interval. 

The control amplitudes at optimization iteration n, 
tii(n), are updated according to 



u{(n) = u 3 k {n-l) + , 



<5$ 



8vJ k {n — 1) 



(15) 



where e is an adjustable scaling factor for the gradients to 
prevent hindering algorithm convergence by making large 
deviations in the control parameters between iterations. 
The algorithm proceeds as follows: 

1. An initial set of control parameters is defined: 

Klf=i(o) 

2. The piecewise constant propagators for each time- 
step are calculated: Uj 

3. The value of the performance functional is calcu- 
lated: $ 

4. If $ > $targ the algorithm exits, otherwise the gra- 
dients for each time-step and control are calculated: 

5. A line-search is used to optimize the step size, e, in 
the gradient direction 

6. The control parameters are updated according to 



7. Loop to Step 2 

To include ringdown suppression and account for distor- 
tions of the ideal control parameters due to the limited 
bandwidth of a high-Q resonator we make several modi- 
fications to this base algorithm. 



B. Optimizing Bandwidth-Limited Controls 

We first resample the control parameters to ensure the 
calculated propagators and gradients accurately reflect 
the distortions of the ideal control fields during periods 
when the applied control voltages are constant. For clar- 
ity, we define a control period as the time At during which 
the applied control voltage is constant, and a evolution 



period as the resampled time At/n s during which the 
control field seen by the quantum system is approximated 
as constant. Concretely, we begin with a set of undis- 
torted controls, {u k }, where j — 1,...,N and resample 
to obtain a new set of undistorted controls {u™}, where 
m = 1, M and M = n s N + n r for n s samples per N 
control periods and n r values of zero appended to the 
waveform to account for pulse ringdown. We then define 
a mapping between the resampled undistorted controls 
and a new set of distorted controls, given by 

the discrete convolution of the undistorted controls with 
the resonator transfer function (JTJ : 



= h 



k*u k 



= E K 



m-l+l~l 



(16) 



Km 



In contrast to the infinite-bandwidth case, the value of 
the control parameters during the m th evolution period 
depends to a certain degree on the value of all previous 
controls. As a result, the gradient of the performance 
functional when control period j is perturbed depends 
on the change in the unitary propagator for all evolution 
periods m > (J — l)n s . As derived in the Appendix, the 
gradients for bandwidth-limited controls are 



M 



5& s—^. 5$ 

fc m>(j — l)n s ^ 



(17) 



where is the convolution of a top-hat function, S, 

with the resonator response: 



(18) 



Km 



The top-hat function, which is formally defined in the 
Appendix, accounts for the resampling of the control 
periods into evolution periods. The convolution of this 
function with the resonator impulse response may be in- 
terpreted as a weighting function of the gradients of all 
evolution periods affected by the perturbation of a given 
control period. The individual gradients of the evolution 
periods are 



5$ 



= -2Re 



~ . At ~ 
P m \i — HkX„ 

n. 



A m Pm 



(19) 



where X m and P m are defined in the Appendix. Note 
that pulse updating is done only for control periods, with 
evolution periods serving as a calculational tool. 

Optimization of the compensation pulse for elimination 
of ringdown is performed as a sub-routine and is not con- 
sidered in the calculation of the gradient direction. How- 
ever, the compensation pulse is taken into account while 
calculating the total pulse fidelity and while performing a 
line search to optimize the step size in the gradient direc- 
tion, e. In practice, the line search is performed by choos- 
ing three values of e, optimizing the compensation pulse 
for each value, evaluating the fidelities of the resulting 
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pulses, fitting to a quadratic, and taking the maximum 
value. The routine for the compensation pulse optimiza- 
tion may be implemented as a two-dimensional search 
over the length and amplitude of an appended final con- 
trol period that minimizes the ringdown. An implemen- 
tation of the GRAPE algorithm for bandwidth-limited 
control is described in Figure [51 



IV. OPTIMIZATION RESULTS AND 
EXPERIMENTAL VERIFICATION 

We performed experiments using a standard irradiated 
fused-quartz sample contained in a Varian E-231 rectan- 
gular cavity on a home-built X-band pulsed ESR spec- 
trometer. The loaded quality factor and impulse response 
function of the cavity were measured by fitting an expo- 
nential to the rising and falling edges of a square pulse 
digitized with a pick-up coil inserted into a cavity iris 
opposite to the sample. Care was taken to only weakly 
couple the pick-up coil to the cavity fields in order to 
disturb the cavity mode structure as little as possible. 

We measured a loaded Q of 8,486 for the Varian cavity, 
giving a ringdown time-constant of r r = 142 ns at a res- 
onance frequency of u!q/2tt = 9.5236 GHz. The r 2 value 
of the exponential fit was 0.9804, verifying that model- 
ing the transient behavior of our high-Q resonator as an 
exponential rise and fall of the field strength is a good ap- 
proximation. We found that a spectrometer deadtime of 
1.2 /is was required to allow the pulse ringdown to decay 
to a value below the spectrometer noise floor. The inho- 
mogeneous phase coherence relaxation time of the sam- 
ple, measured by observing a spin echo, was T 2 * ~ 250 ns, 
preventing us from observing a useful FID of the quartz 
sample when using an unoptimized pulse. 

To observe the quartz FID we optimized a bandwidth- 
limited OCT pulse that performs a high-fidelity ir/2 ro- 
tation about the x-axis robust to variations in the static 
and microwave fields. The resonator model discussed in 
Section [TT] was included in the optimization, as well as 
ringdown suppression. The Hamiltonian used for opti- 
mization, in a frame rotating at the nominal electron 
Larmor frequency, was 

H(Au,ui) = -Auht z + -wiA(t)cT x , (20) 

where Aw is a resonance offset parameter representing 
static field inhomogeneity in units of rad/s, w\ is a scal- 
ing factor of the nominal Rabi frequency, wi ino mj repre- 
senting microwave field inhomogeneity, and A{t) is the 
time-dependent amplitude modulation representing the 
OCT pulse, with range {— u)\ tXlola , Wi,nom}- The spin dy- 
namics was calculated by taking a convex operator sum 
over a uniform classical probability distribution of field 
inhomogeneities, P(Aw,a;i) 37]. 

The pulse was defined piecewise constant over 100 
steps of 10 ns each, giving a total length of 1 fis. The 



nominal Rabi frequency, uji. nonl /2ir = 5.26 MHz, was de- 
termined by Fourier transforming the result of a standard 
spin-echo Rabi oscillation experiment at a microwave 
power of 4 Watts and identifying the dominant frequency. 
The pulse was optimized over a microwave inhomogene- 
ity of bJi = (0.95, 1, 1.05) and a static field inhomogene- 
ity of Aa;/2vr = {-2 MHz, 2 MHz} in steps of 250 kHz. 
The pulse optimization took roughly ten minutes on a 
standard laptop computer and resulted in an average fi- 
delity over the distribution of $ = 0.9905. The resulting 
pulse profile and response over an extended distribution 
is shown in Figure [3] The free-induction decay (FID) of 
irradiated fused-quartz resulting from application of the 
pulse in the Varian cavity is shown in Figure^ along with 
the digitized pulse profile measured through the pick-up 
coil. A spectrometer deadtime of 75 ns was required to 
allow the small oscillations shown at the end of the digi- 
tized pulse to decay. This deadtime was included in the 
pulse optimization as an additional control period of zero 
amplitude, allowing observation of the full FID at high Q 
and verifying the simulated behavior of the bandwidth- 
limited OCT pulse 



V. CONTROLLABILITY WITH LIMITED 
BANDWIDTH 

For accurate spectroscopy and high-fidelity quantum 
information processing, control sequences must drive 
transitions and excite coherence over a range of frequen- 
cies given by the coupling structure and any uncertainties 
of the system Hamiltonian. A common solution, based 
on linear response theory (39j . is to require the Fourier 
spectrum of the control pulse to contain significant con- 
tributions from all frequencies present in the quantum 
system [4(], Sl|. The principle, and limitations, of this 
approximate solution may be succinctly demonstrated 
using a description of the pulse trajectories in fc-space 
(42| . In this description, each frequency u> present in the 
drift Hamiltonian of the quantum system has an associ- 
ated wavenumber, k = ut, whose value is modulated by 
the control Hamiltonian during the length of the applied 
pulse. 

The accuracy of predicting spin response through 
Fourier analysis depends primarily on the validity of two 
approximations. The first is that any rotations that oc- 
cur during a period where the pulse waveform is constant 
are small, such that the full spin response given by the 
exponential of the Hamiltonian generating the motion 
may be truncated to first order. The second is that the 
amplitude of the control Hamiltonian during any pulse 
period be significantly greater than the amplitude of the 
drift Hamiltonian, such that the axis of rotation may 
be taken as being completely determined by the control 
Hamiltonian. When these two criteria are satisfied, exci- 
tations of the quantum system as a function of frequency 
may be written simply as the Fourier transform of the 
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time-dependent control amplitudes, x(t): 



S{oj) =-i x{t)e- lk dt, 
Jo 



(21) 



where t p is the pulse length and S(lu) is the observable 
spin response generated by the controls, which are taken 
for simplicity to be amplitude modulated only. If we now 
write the controls as a sum of Fourier components, 



(22) 



we see that the only contributions to the observable sig- 
nal are from Fourier components, n, where i„ such 
that the corresponding vector in fc-space may be effec- 
tively refocused during the pulse. Within the linear re- 
sponse approximation, then, pulse design is frustrated by 
the narrow-band filtering property of a high-Q resonator. 
However, the ability to numerically optimize pulses al- 
lows us to move into a regime where linear response is no 
longer valid. 

For general rotations there is significant mixing of the 
various fc-vectors that leads to a complex spin response. 
Also, a general axis of rotation is not given entirely by the 
control Hamiltonian, but by the vector sum of the con- 
trol and drift Hamiltonians. These effects are often small, 
but may be used to generate spin response deemed inac- 
cessible by linear response. By solving the equations of 
motion exactly under an accurate system model, without 
making approximations, we can find bandwidth-limited 
pulses that retain the same degree of controllability as 
for infinite bandwidth pulses, albeit with a reduction in 
efficiency due to the inability to directly address all tran- 
sitions in the system. 

We first consider an ensemble of uncoupled spins with 
Larmor frequency, luq, off-resonance an amount Auj = 
luq — uj t from a control field applied at uj t . The Hamilto- 
nian in a frame rotating at uj t is given by eq. (|20[) , where 
the term proportional to a z is the drift Hamiltonian, Hd, 
and the term proportional to a x is the control Hamilto- 
nian, H c . Generating a computationally universal set of 
quantum operations requires the ability to generate all 
elements of the Lie algebra spanning the Hilbert space 
43]. For uncoupled spins, the Lie algebra is SU(2), with 
basis operators {I,a x ,a y ,a z }. The control algebra gen- 
erated by a given set of Hamiltonians may be computed 
by taking successive Lie brackets to all orders 
For Hd and H c considered here [Hd, H c ] cx 
the generated control algebra is identical to SU(2), in- 
dicating universal control of the system by appropriate 
application of the given Hamiltonians. 

The efficiency of control is difficult to quantify in gen- 
eral, but may be posed as the non-commutativity of ef- 
fective Hamiltonians that may be generated during the 
pulse between incremental time periods. For SU(2) we 
may take a simple geometric view, posing the problem 
of efficiency as the maximum angle of rotation that may 



imputed 

mm. 

such that 



be generated over an incremental time period or, equiv- 
alently, the maximum angle between effective Hamilto- 
nian vectors that may be generated from one instant in 
time to the next. Assuming an exponential model of the 
pulse transients for simplicity, and taking A(to) — and 
uj\ = 1, we examine the non-commutativity of 



H cS (t ) = -Au}a z , 



(23) 



and 



H cS (t Q + 6t) = -Auja z + -A{t + 5t)(l-e-^- 5t )a x . (24) 

The controllability of the system has not changed, pro- 
vided Q is finite, as these Hamiltonians generate the same 
Lie algebra as for infinite bandwidth. The angle between 
the two Hamiltonian vectors, 0, is given by 



= tan 



A(t Q + St) 
Auj 



(1 



(25) 



The size of the generated incremental rotation depends 
on the ratio of the target control amplitude to the reso- 
nance offset and the value of Q, with increasing efficiency 
for larger control amplitudes and smaller Q. Given that 
control amplitude is proportional to \fQ, there is an in- 
herent trade-off between sensitivity and the efficiency of 
generating rotations that must be considered. However, 
given sufficient time and a sufficiently accurate system 
model, pulses may be optimized that address transitions 
that are arbitrarily far off-resonance without containing 
Fourier components anywhere near the transition fre- 
quency. In the following section we provide examples, 
in the context of anisotropic hyperfine coupled electron- 
nuclear spin systems, of pulses that significantly exceed 
the limitations imposed by linear response theory. 



VI. APPLICATION OF BANDWIDTH-LIMITED 
CONTROL TO ELECTRON-NUCLEAR SPIN 
SYSTEMS 

We consider a one electron, one nuclear (le-ln) spin 
system in a strong external magnetic field Bq = Bqz 
with an anisotropic hyperfine interaction coupling the 
two spins. The drift Hamiltonian is given by 



H d 



1 



1 



, , C _L- I , ,71. _f. 



-w zz a e z a z 



1 



(26) 



where ui e z and uj z are the strengths of the electron and 
nuclear Zeeman interactions, respectively, with the static 
field, and uj zz and ui zx are isotropic and anisotropic com- 
ponents of the hyperfine interaction. Examples of such 
systems include solid-state organic molecules with a lo- 
calized free radical 0, [47} , Phosphorus donors in silicon 
(48j , and nitrogen- vacancy defect centers in diamond [4!| . 

The electron Zeeman interaction is dominant and de- 
fines the principle coordinate system of spin quantiza- 
tion. In terms of this coordinate system, the eigenstates 
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of the nuclear spins are given by the vector sum of the 
nuclear Zeeman interaction and the spin dependent local 
field generated by the anisotropic hyperfine coupling to 
the electron: 



|1) = |ta o )=sin0 t |tt)+cos0t|t+>, 
|2> = |tai)=cos0 t |tt}-sin0 t |t+), 
|3)H+/3i}=cos0;|+t)-sin0;|W), 
|4> = |+A)}=sin^|+t>+cos^|W), 



(27) 



where and #4. determine the non-commutativity of the 
resulting eigenstates and are given by: 



6^ = tan 



= tan 



Uzz - U' z 

~^zx 
UJzz + u) 7 ; 



(28) 



Due to the non-commutativity of the local fields seen 
by the nuclear spin when the electron spin is in the spin- 
up versus spin-down state, universal control over the en- 
tire spin system (generalizing to le-Nn systems) may be 
achieved by implementing only the electron o e x genera- 
tor and allowing free evolution under Hd [5Cj | . This may 
also be seen by noting that the transition probability be- 
tween any pair of non-degenerate states associated with 
the 0% operation is non-zero [5l|. We may thus use the 
nuclear spins as a quantum processor, with the electron 
spin acting as an actuator element to allow for fast quan- 
tum operations on the processor and providing a means 
for efficiently transferring information in a node-based 
quantum information processor design [12, [5(1 [HJ ■ 

To demonstrate that universal control via electron-only 
modulation may still be achieved when using a high-Q 
resonator we numerically optimized a pulse which imple- 
ments an electron spin flip, 



U d = e~^' 



(29) 



for a le-ln sample in a resonator with bandwidth smaller 
than the size of the hyperfine interaction. The ability 
to achieve such an operation implies full controllability 
of the spin system with electron-only modulation. The 
control Hamiltonian was taken to be a time-dependent 
amplitude modulation of the electron spin only, 



1 



H c {t) = -A(t)oZ®r 



(30) 



The pulse was optimized using the algorithm outlined 
in IIII Bl using a model of the resonator given in [XT] with 
a Q of 10,000 (BW w 1.2 MHz). The drift Hamilto- 
nian parameters were taken from [H3|: uj z /2tt = 11.885 
GHz, w"/2tt = 18.1 MHz, uj zz /2tt = -42.7 MHz, and 
lj zx /2tt — 14.2 MHz. The carrier frequency of the pulses 
was set resonant with the 1-4 transition, lu /2tt = 11.909 
GHz. The nominal Rabi frequency was taken to be 
^nom/27r = 100 MHz. The resulting pulse is shown in 



Figure [5a] A pulse time of 5 [is was chosen for conve- 
nience to provide sufficient time to easily achieve the de- 
sired operation. We expect solutions for shorter pulse 
times to exist, but have not systematically addressed 
the minimum time needed to achieve the operation for 
a given set of parameters. The final simulated average 
gate fidelity of the pulse was 0.9901 even when, as shown 
in Figure [5bl all significant frequency modulation is much 
less than the hyperfine splitting. 



VII. CONCLUSIONS 



Integrating the impulse response function of a high-Q 
resonator into an optimal control theory pulse design al- 
gorithm allows higher sensitivity and signal-to-noise ratio 
in measurements without sacrificing the ability to per- 
form a universal set of quantum operations with high 
fidelity. We focused on a particular application in pulsed 
electron spin resonance where universal control of the 
quantum state of a nuclear spin may be achieved by 
microwave-only modulation of an anisotropic hyperfine 
coupling to an electron spin. We demonstrated that the 
limits imposed by linear response theory may be vastly 
exceeded when using numerically optimized pulses, al- 
lowing universal control of the joint electron-nuclear spin 
system even when the hyperfine coupling strength is sig- 
nificantly greater than the resonator bandwidth. These 
solutions rely on the non-commutativity of the effective 
Hamiltonians generated throughout the pulse and the 
ability to accurately engineer the exact response of a 
quantum system to complex modulations. Experimen- 
tal data verified the accuracy of our system model and 
the practicality of our presented optimization method. 

We expect the methods presented in this work to find 
application in a broad range of fields. In addition to 
spectroscopic and imaging applications in chemistry, bi- 
ology, and materials science for samples with a small 
number of spins, we expect the application of high-Q 
control techniques to molecular thin film samples to be 
particularly compelling due to their relevance to building 
a large-scale multinode quantum information processor. 
The described techniques may also find application in hy- 
brid quantum systems that aim to use spin ensembles as 
memory elements for microwave photons (53l . l54j . 

We note that the optimized pulses work well only for 
the specific set of system parameters defined at the time 
of optimization and the performance is sensitive to any 
change or mis-setting of the system parameters when run- 
ning an experiment. The pulses may be made robust to 
variations in system parameters, but with a correspond- 
ing increase in pulse time and complexity [37], HH, IBtij . 
The efficiency of control when both robustness and high- 
fidelity are required of pulses to be used with a high-Q 
resonator requires further investigation. 
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VIII. APPENDIX: DERIVATION OF 
DISTORTED GRADIENTS 

We begin with a set of undistorted controls, 
where j = 1,...,N and rcsample to obtain a set of 
undistorted controls {u m }, where m = 1, ...,M and 
AI = n s N + n r for n s samples per N control periods 
and n r values of zero appended to the waveform to ac- 
count for pulse ringdown. In this derivation, we set k = 1 
for clarity. The distorted controls, {v m }, are then given 
by the discrete convolution 



The average gate fidelity may then be written as 



$ = (U d \U M ...Ui)(Ui...U M \U d 



(31) 



(32) 



where the m th propagator is calculated from the distorted 
controls as 



Un 



We are interested in calculating the gradient of the fi- 
delity with respect to a perturbation of the undistorted 
controls {u J }, given by the usual product and chain rules 
of derivation as 



m V In \t~i SUm 

m— 1 \ 



U d \U M ...Ui) (u x 



.U x (U 1 ...U M \U d ) + 



SUrn 



(34) 



■U M \U d 



By defining a forward propagator, X m = U m ...U\, a 
backward propagator, P m = U m+l ...U' M U d , and noting 
that only propagators with m > (J—l)n s will be modified 
by a change in the j th control period, we may simplify 
this expression to 




X m -i\P n 



(35) 



In the following, due to resampling the control period, 
we define a top-hat function, n s ), to account for the 
perturbation being constant during the entire j th control 
period 

n s ) = 1 for (j — l)n s < m < jn s + 1 (36) 



n s ) = otherwise 



(37) 



The derivative SUm/Su^ may then be calculated by defi- 
nition as 

5U m .. U m (v m +E{j,n s )5ui)-U m (v m ) ,„ Q . 
-j — - = lim — : (38) 

We apply a small perturbation to the j th control, vP — > 
u J + Su^ and determine the resulting m th propagator: 



t/ m (F™ + ~(j, n 8 )6u*) 



3Ci,n,)<yu^)Hi 



(39) 



Assuming At to be small, we may approximate H d and 
Hi as commuting, allowing us to keep only the first order 
term of a BCH expansion of the perturbed propagator 



C/ m (F" + S(i,n s K) 



1 S(j,n B )Su J H t 



(40) 



If we now assume the perturbation Su^ to also be small, 
we may make a first order series approximation of the 
exponential, to give 

U m (v m + Z(j,n s )uJ) = 
(33) U m (v m )(I - i—Hjj m h m - l+1 E(j, n s )8ui + 0(At 2 ) 



Ki 



(41) 



We may now plug this expression back into (|38|) to obtain 



SUm 

8vj 



At 



— R x U m y2h m - l+1 E(j 7 n s ) (42) 



Km 



which, when plugged into (|35j) gives 

M 



5$ 



Sv r 



(43) 



m>(j — l)n s 

where £, m (j) is the convolution of the top-hat function 



rtf) = X> m-,+ls & n «) 



(44) 



Km 



This function may be interpreted as a weighting of the 
change of the fidelity due to a small change Sv m of the 
pulse parameter at the m th evolution time step induced 
by the perturbation Su J 



6$ 

5v m 



= -2Re 



Pm\i H\X m ) (X m \Pn 



(45) 



Note that when the controls are undistorted, the gradient 
weighting function is equivalent to the top-hat function, 
giving 



S<S> 



E 



5$ 



(46) 



»>(j'-l)n.+l 



which for short time steps and small perturbations, re- 
duces to the undistorted gradient derived in [33| 

H = -2Re [(PjliAmXj) {X 3 \P 3 )\ (47) 
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FIG. 1: Transient behavior of a resonator circuit. A general resonant transmission circuit may be modeled as 
a series tuned RLC circuit capacitively coupled to a time- varying voltage source [22| (inset). The quality factor (Q) 
of the resonator is given mainly by the circuit resonance frequency, luq = 1/^/LCt, the coil inductance,!/, and the 
coil resistance, r: Q ~ uj^L/r. The circuit impedance is matched to Rq — 50 Ohms by varying the capacitances Ct 
and Cm- For high-Q resonators, the dominant transient response of the resonator to a square pulse input of voltage, 
Vg, (bold line) is an exponential rise (dashed line) and subsequent ring-down (dotted line) of the coil current, it,, 
and resulting magnetic field. The ringdown may be suppressed by application of a phase inverted compensation 
pulse at the end of the square pulse to drive the coil current to zero. The characteristic transient ringdown time 

without compensation, r r = Q/u>o, is denoted by a cross. 
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FIG. 2: Schematic of optimization algorithm. The bandwidth-limited GRAPE pulse optimization algorithm 
proceeds in a similar manner as for undistorted controls (33| , with the resonator transfer function included in the 
calculation of the pulse propagator, average gate fidelity, and gradients. The notation is explained in the main text, 

with k — 1 for simplicity. An initial guess of control amplitudes is resampled and convolved with the resonator 
transfer function to yield a distorted set of control amplitudes including ringdown. A compensation pulse period is 
then optimized in a sub-routine and appended to the waveform to yield a distorted set of control amplitudes with 
minimized ringdown. The propagator for each evolution period of the distorted control amplitudes is then calculated 
and the average gate fidelity computed. The set of propagators is used to calculate the gradient direction of the 
fidelity with respect to the undistorted controls. A line search is then performed to optimize the step-size in the 
gradient direction. At each step of the line-search the sub-routines to calculate the updated distorted control 
amplitudes and corresponding compensation pulse are called, accounting for the compensation pulse not being 
included in the gradient calculation. The undistorted controls are then updated according to the gradient direction 
and step-size and the process iterates until a desired value of the fidelity is achieved. 
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FIG. 3: Bandwidth-limited OCT pulse robust to static and microwave field inhomogeneities. (a) 

Time-domain profile of a transients optimized OCT pulse with ringdown suppression implementing a ir/2) x rotation 
robust to variations in the static and microwave field strengths. The dashed line represents the undistorted controls 
and the solid line represents the control fields seen by the spin system after transmission through a resonator with Q 

= 8,486. The control amplitudes are normalized to a nominal Rabi frequency of wi inom /27r = 5.26 MHz. 
Optimization parameters and further details are discussed in the main text. Note that this pulse is phase- refocused, 
in that all spins in the sample are rotated with the same phase, (b) Simulated pulse fidelity over an extended range 
of static and microwave field inhomogeneities, demonstrating the robustness of the pulse. 
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FIG. 4: Digitized pulse and free-induction decay (FID) of irradiated quartz. The transients optimized 
OCT pulse shown in the previous figure was applied to an inhomogeneously broadened solid-state sample of 
irradiated fused-quartz in a rectangular cavity with Q = 8,486 at a resonance frequency of 9.5236 GHz. The black 
solid line for t < shows the pulse profile digitized through a pick-up coil inserted in the cavity. The digitized profile 
closely matches the calculated profile with the only ringdown being a small oscillation which decays after roughly 75 

ns. The resulting quartz FID is shown as a blue solid line for t > and was acquired after a 75 ns spectrometer 
deadtime. The static field was moved roughly 2 MHz off-resonance (still within the high-fidelity operation regime of 
the pulse, as shown in Figure [3b]) to emphasize the shape of the FID. For comparison, a digitized square pulse of 
length 1 [is is shown as the black dashed line. In a separate measurement, the necessary spectrometer deadtime in 
the absence of ringdown suppression was determined to be 1.2 /is (shown as a dotted vertical line), which would 
prevent the detection of a significant portion of the FID. The separate plots have been scaled for visual clarity. 
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FIG. 5: Optimized le-ln pulse at high-Q. (a) The undistorted (dashed line) and distorted (solid line) control 
amplitudes of a pulse designed to perform a tt rotation of the electron spin in a resonator with Q = 10,000. Access 
to such an operation guarantees universal control of the nuclear spin via electron-only modulation [50( . The pulse 
consists of 500 time steps of At= 10 ns each for a total length of 5 /zs, with wi inom /27r = 100 MHz. The simulated 
average gate fidelity is $ = 0.9901. (b) The single-sided amplitude spectrum of the undistorted pulse (dotted line) 
and the distorted pulse (solid line) filtered by the resonator admittance function (bold solid line). The transitions 
necessary to achieve the desired operation were separated by an amount {\u>23 — ^14 1 = 51 MHz) much greater than 
the bandwidth of the resonator (w 1.2 MHz), demonstrating that linear response may be greatly exceeded by pulses 
of high complexity optimized under a sufficiently accurate system model. 



